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Abstract 



I We take into account the sigma meson cloud effect in the meson cloud model to 

\^ ' calculate the distributions of light flavor sea quarks in the proton. Our calculation 

On ■ gives a better description of the data for d{x)/u{x). We also provide a picture that 

the probability of finding a physical proton in a Fock state |A^a;) is reasonable small 
O '• with a smaller cutoff Aa,. 
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^ ' By now, it has been experimentally established that the light-flavor sea 

quarks d{x) and u{x) in the proton differ substantially [1,2,3,4]. The violation 
of the Gottfried sum rule, first found by NMC [1], indicated that D = jQ[d{x) — 
u{x)] dx = 0.148 ± 0.039, and similar result for d{x) — u{x) was obtained by 
HERMES [2]. In Drell-Yan experiments, d/u was determined to be more than 
2 at X = 0.18 by NA51 [3] and the Bjorken-x dependence (moment fraction 
dependence) of the ratio has been measured by E866/NuSea [4]. 

While no known symmetry requires d(x) equal to u{x), the large asym- 
metry was out of naive expectations. The sea of light flavor quark- ant iquark 
pairs produced perturbatively from gluon splitting is flavor symmetric. Thus 
the large asymmetry requires a nonperturbative origin. It was flrst predicted 
by Thomas [5], and later suggested by Henley and Miller [6], and Kumano 
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and Londergan [7] , that including the effects of virtual mesons (recognized as 
meson cloud model) can naturally account for the asymmetry. In previous ar- 
ticles applying meson cloud model to explain the light flavor sea asymmetry in 
the nucleon, a good description of the data was obtained for the distribution 
functions [8]. The dominant role is played by the pion, which provides that 
d{x) /u{x) either increases monotonically with x or turns back towards unity 
too slowly [8,9]. Many explanations have been applied to attack this prob- 
lem [10,11,12,13], such as effects of A [10,11,12], the influence of the Pauli 
exclusion principle [12], adjustment of parameters [13], but none of these ex- 
planations provides a satisfactory description of the ratio d{x)/u{x). However, 
Alberg, Henley and Miller [14] found that the inclusion of the effect of the 
isoscalar vector meson uj, with a coupling constant 5'5/47r ^ 8.1, allows a good 
description of the data. As mentioned in that article, adding the a effect will 
likely improve the description of the data. Here, we follow this suggestion and 
take into account the sigma meson effect in the meson cloud model. 

The existence of the a meson has been obscure for many years, while many 
theorists and experimentalists were searching for a which may play an im- 
portant role in nuclear physics, because the a meson can provide reasonable 
middle-range nuclear force [15]. Sigma meson was flrst introduced as the chi- 
ral partner of pion in the linear representation of chiral symmetry, which is 
an important ingredient in modern hadron physics. It also plays an impor- 
tant role in spontaneous chiral symmetry breaking to understand the present 
spectroscopy of hadrons. However, it was not well established experimentally 
mainly due to the negative results of extensive analysis of the vr-vr scatter- 
ing phase shift. While recent reanalysis of the vr-vr scattering phase shift [16] 
strongly suggested a pole of the a-meson. Though we do not have a conflrmed 
conclusion that the sigma meson exists as a real physical particle, it is suffi- 
cient for considering its effect to calculate d{x) — u{x) and d{x) /u{x) in the 
proton in the meson cloud model. 

In the meson cloud model [7,8,9], the nucleon can be viewed as a bare nu- 
cleon (core) plus a series of baryon-meson Fock states which result from the 
fluctuation of nucleon to baryon plus meson N — *■ BM (a bare baryon sur- 
rounded by a meson cloud). The model assumes that the life-time of a virtual 
baryon-meson Fock state is much longer than the interaction time in deep in- 
elastic scattering (DIS) or Drell-Yan processes, thus the quark and anti-quark 
in the virtual baryon-meson Fock states can contribute to the parton distri- 
bution of the nucleon. These non-perturbation contribution can be calculated 
in a convolution between the fluctuation function, which describes the micro- 
scopic process — > BM, and the quark (anti-quark) distribution of hadrons 
in the Fock states \BM). Here we provide the usual formula. First the physical 
proton wave function is composed of the following Fock states 



2 



|P) = |P),_ + [ dyd'k^^BMiy, kl) \B{y, k^), M(l - y, -k^)) . 

BM •' 

Here Z is the wave function renormalization constant, '^BM{y, k^) is the prob- 
abihty amphtude for finding a physical nucleon in a state consisting of a baryon 
B with longitudinal moment fraction y and transverse momentum k±, and a 
meson M with longitudinal moment fraction 1 — y and transverse momentum 
— k±. Therefore, the quark distribution functions q{x) in the proton are given 
by 

q{x) = Zq^,,,{x) + 5q{x) , (1) 

with 



H^) = Y.[j-fMBiy)qMi-) + j -fBM{y)qB{-)] , (2) 

mb{ y y i y y 

and 



fMBiy) = fsMil - y) ■ (3) 

As shown in [11], the wave function renormalization constant reads 

1 1 
^ = [1 + E / fBM{y)dy]~' ^ 1 - E / fBM{y)dy . (4) 

BM Q BM Q 

The number of each type of meson, um, is obtained by integrating fMsiy) 
over y. Then we have Z = 1 — Y1,m^m- The splitting function is 

oo 

fBM{y) = j |vl>BA/(y,kl)pd2k^. (5) 


We use time-ordered perturbation theory (TOPT) in the infinite momentum 
frame (IMF) to calculate this function, which was given by [8] as 

r I ^_ 1 mNimB \GM{y,'k\)?\VmF? 

where 
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y 1-2/ 



(7) 



is the invariant mass of the final state, and an exponential form for the cutoff 
is [11,15] 

GMiy, = exp( —2 ) , (8) 

which is required to respect Eq.(3). Here Am is a cutoff parameter for each 
meson. 

The vertex function ViMpd/ik^) depends on the effective interaction La- 
grangian that describes the fluctuation process — BM. Here we use the me- 
son exchange model for hadron production [15] and the linear sigma model [17] 
to perform the calculation. 

In a full calculation, we should include all kinds of mesons and baryons. 
While the probability of baryon-meson fluctuation should decrease with the 
invariant mass of the baryon-meson Fock state increasing, we can neglect the 
effects of Fock states with higher invariant mass. We just include speciflcally 
71, a and u mesons with the nucleon here. It has been discussed that the effect 
of the p meson and the intermediate A can also contribute a lot each. But 
the p meson increases d{x) — u{x), whereas the intermediate A decreases it, 
so these effects tend to cancel each other. And we also omit the effect of t] for 
its small coupling constant. Thus it is safe to focus on the effects of tt, a and 
u meson clouds here. 

First, we focus on the sigma meson. From the Bonn meson-exchange model 
for the nucleon-nucleon interaction [15], the effective interaction Lagrangian 
for the scalar meson a is L = gipaip. Thus we can get the vertex function for 
a to obtain the fluctuation function as 

f^^(y) = f-.-ii^, /dkiia(y,ki)P . (9) 

87r2y2(i_^)j [m^-M^^(?/,k^)]2 

Now, we should know the coupling constant g^, the mass and the cutoff 
A.. 

As mentioned above, sigma meson is introduced to the hadron physics as 
the chiral partner from the view of chiral symmetry. In the SU{2) linear a 
model [17], the effective Lagrangian is 



4 



2.0 




1.6 - 



A_ = 1.0 

a" = 0.88 

a] = 0.83 



■ E866/NuSea 
• HERMES 



i 



0.0 



0.1 



0.2 
X 



0.3 



0.4 



0.5 



Fig. 1. Comparison of our meson cloud model with data [2,4] for d{x) — u[x). The 
thick solid curve is for = 0.88 GeV, which gives a best description here. 



Even though we do not have direct information about Qa, = dn is imposed 
by the Unear a model According to all kinds of calculations to fit different 
experimental data [15,18,19], the coupling constant gl/4:H = gl/4:7i = 13.6 is 
taken in our numerical computations. In 2002 Particle Data Group (PDG), 
/o(600) or a appears below 1 GeV mass region. Recently it was found that 
m„ = 585 ±20 MeV by a reanalysis of tttt scattering phase shift [16]. Similarly, 
by analysing tttt production processes [16], we obtain mo- = 580 ± 30 MeV for 
pp central collision. Thus, we will set the mass as 600 MeV in the following 
calculations. It is generally believed that the cutoff A^. value is in a range 
around 1 GeV as a phenomenological parameter. Here we will examine the 
effect of varying Ao- value in the range 1.0 < Ao- < 1.3 GeV. 

We have made some detailed discussions about the a meson cloud above. 
The vr meson cloud effect has been discussed enough in previous articles 
[6,8,14,20], here we just set the parameters as = 139 MeV and A^ = 
(0.88 ± 0.05) GeV, which are chosen to reproduce the integrated asymmetry 
D = /o^[(i(x) — u{x)] dx [1,2] as shown in Fig. 1. 

Moreover, we should discuss the functions of qM{.x) and qsix). Those for 
the nucleon and pion are measured, but the quark distribution functions of 
the a meson and u meson are unknown. It has been traditional [21] to assume 
that the structure function of the p and n are the same. It is suggested by the 



L = gip{a + -f5T ■ 7i)ip . 



(10) 
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Fig. 2. Comparison of a harder bare nucleon sea quarks (thick soUd Hne) with a 
traditional bare nucleon sea quarks (dashed line). The thin solid line is only the vr 
contribution to the ratio d{x)/u{x). 

bag model [14] that the structure functions of the cj, p, and tt are the same. 
Due to that the a meson operator is adopted as -^(uu + dd), there is some 
reason to assume that the quark distribution functions of the a meson is equal 
to that of the n. Thus we use [22] 

xg_(a;)= 0.2(1 (11) 

for the valence and sea quark distribution functions of all mesons considered 
here. 

The bare nucleon sea is traditionally parametrized as 



xq,^Jx) = 0.11(1 -xY'-^ 

which is recognized as Holtmann's parametrization of the bare nucleon sym- 
metric sea. While Alberg and Henley [20] used a harder distribution for the 
bare sea quarks of the form found in the determination of the gluon distribu- 
tion [24] 

0.0124a;-°-^^(l -a;)^-^ (13) 
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Fig. 3. Comparison of our meson cloud model with data [4] for d{x)/u{x). The 
thick solid curve shows our result only considering the vr meson contribution. The 
shadowed area shows the effect of adding a meson and varying Ao- value in the range 
1.0 < A<^ < l.SGeV. 

Their calculations indicate that the harder distribution for the bare sea quarks 
gives a better description of d{x)/u{x) as shown in Fig. 2. In the following 
calculations, we will use only the harder bare nucleon sea. Also we can see that 
the change is not enough to explain the data, which guarantees the necessary 
to consider other flavor symmetric contribution. The inclusion of isoscalar 
meson uo with a reasonable coupling constant produces a similar improvement 
in agreement between theory and experiment [14], so does the inclusion of 
isoscalar meson a. 

The results of our calculations for light flavor sea quarks are shown in Fig. 3 
and Fig. 4. In Fig. 3, we examine the effect of varying A^- value in the range 
1.0 < Act < 1.3 GeV. The thin solid curve on the upper of the shadowed area 
is the result of setting Act = 1.0 GeV, while the bottom curve corresponds 
to Act = 1.3 GeV. It is clear that adding sigma meson effect here gives really 
better description of the data. Our calculations illustrate that the larger value 
of Act tends to give small values of the ratio d/u and that decreasing Act causes 
the maximum value of d/u to be larger and to appear at a higher value of x. 

In Fig. 4, we illustrate the effects from three mesons with different cut- 
off values. As calculated in [14], we set qIj/^t^ as 8.1, which is from fit- 
ting to dispersion relation descriptions of forward nucleon-nucleon scattering 
[15,19,25]. The shadowed area shows the results of varying A^ value in the 
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Fig. 4. Comparison of our meson cloud model with data [4] for d{x)/u{x). The thick 
sohd curve shows our result if contributions from the a and lo clouds are omitted. 
The dashed curve shows adding lo cloud contribution with = 1.0 GeV. The 
shadowed area shows the contributions of three mesons, and the effect of varying 
value in the range 1.0 < < l.SGeV. 

range 1.0 < < 1.5 GeV with at 1.0 GeV on the left and 1.3 GeV on 
the right. The dashed curves omit the a contribution and only consider the 
value of as 1.0 GeV. When the value of A^ is 1.5 GeV, the ratio d{x)/u{x) 
hardly changes with adding a meson cloud, which is consistent with [14], in 
which is set as 1.5 GeV to obtain a good description of data by considering 
only 71 and uj effects. In order to have a deep understanding of this, we present 
the numbers of the related mesons um in the proton with the different cutoff 
values in the two works, which are shown in Table 1. The large value of cutoff 
leads to a large value of meson number in the proton. Apparently, the too 
large resulted from the large value of A^ = 1.5 GeV should not be the 
real picture of the proton. For this reason, A^^ favors the small value of the 
examined range. Adding a cloud effect modifies the parameters describing the 
omega-nucleon interaction. And such a revision provides a picture that the 
probability of finding a physical proton in a Fock state \Nuj) is reasonably 
small. We also realize that the same rule holds true for the a-meson from the 
comparison of the two different values of A^-. 



Table 1. Parameters and meson numbers in the proton 





A^ 


A. 


A,^ 








Our work 


0.88 


1.0 ~ 1.3 


1.0 ~ 1.5 


0.175 


0.023 ~ 0.078 


0.063 ~ 0.671 


Alb erg-Henley's work 


0.83 


no sigma 


1.5 


0.150 


no sigma 


0.671 



In summary, the inclusion of the a meson cloud effect brings a better de- 
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scription for d{x) /u{x) in the proton and also provides a picture of a reasonable 
small in the proton. Thus we conclude that the inclusion of both isoscalar 
meson sigma and omega cloud effects can give an improved description of the 
experimental data for d{x) /u{x) in the meson cloud model. 

This work is partially supported by National Natural Science Foundation 
of China under Grant Numbers 10025523 and 90103007. 
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